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Abstract. If D is a rational polygon, then the associated rational billiard 
table is given by 0(D). Such a billiard table is well understood. If F is a 
closed fractal curve approximated by a sequence of rational polygons, then 
the corresponding fractal billiard table is denoted by Q{F). In this paper, we 
survey many of the results from [LapNiel-3] for the Koch snowflake fractal 
billiard f2(KS) and announce new results on two other fractal billiard tables, 
namely, the T-fractal billiard table (see [LapNie6]) and a self-similar 

Sierpinski carpet billiard table Q{Sa) (see [CheNie]). 

We build a general framework within which to analyze what we call a 
sequence of compatible orbits. Properties of particular sequences of compatible 
orbits are discussed for each prefractal billiard Q(KSn), ^{^n) and Q(Sa,n), 
for n = 0, 1, 2 • • ■ . In each case, we are able to determine a particular limiting 
behavior for an appropriately formulated sequence of compatible orbits. Such 
a limit either constitutes what we call a nontrivial path of a fractal billiard 
table f!(i^) or else a periodic orbit of Q{F) with finite period. In our examples, 
F will be either KS, S or Sa- Several of the results and examples discussed 
in this paper are presented for the first time. 

We then close with a brief discussion of open problems and directions for 
further research in the emerging field of fractal billiards. 
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1. Introduction 



This paper constitutes a survey of a collection of results from [LapNiel, LapNie2, 
LapNieS] as well as the announcement of new results on the T-fractal billiard table 
(see [LapNie6]) and a self-similar Sierpinski carpet billiard table ri(S'a) (see 
[CheNie]). 

In §§2 and 3, we survey the necessary background material for understanding the 
remainder of the article. More speciflcally, in §2, we introduce the notion of a ratio- 
nal polygonal billiard, a flat surface determined from a rational polygonal billiard 
and discuss the consequence of a dynamical equivalence between the billiard flow 
and the geodesic flow.^ This dynamical equivalence allows us to express an orbit of 
a rational billiard tabic as a geodesic on an associated flat surface, and vice-versa, 
with the added benefit of being able to determine the reflection in certain types of 
vertices of a rational billiard table. Furthermore, in §3, we provide additional back- 
ground material from the subject of fractal geometry necessary for understanding 
the construction of the Koch snowflake KS, T-fractal and a Sierpinski carpet 
5a, as well as particular orbits and nontrivial paths. 

We then combine the background material presented in §§2 and 3 to analyze the 
prefractal billiard tables rt{KSn), ^{^n) and ^{Sa,n), for n = 0, 1, 2 • • • . We begin 
by providing a general language for prefractal billiards and subsequently focus on 
determining sufficient conditions for what we are calling a sequence of compatible 
periodic orbits. While §§4.2-4.4 contain specific results and specialized definitions, 
there is an over-arching theme that is more fully developed in §5. 

In addition to providing a general language within which to analyze a fractal 
billiard, we discuss in §§5.2-5.4 how one can determine well-defined orbits of fl{KS), 
il(^) and fl{Sa), as well as nontrivial paths of fl{KS) and Sl{^) that connect two 
elusive points of each respective billiard. Relying on the main result of [Du-CaTy], 
the second author and Joe P. Chen have shown that it is possible to determine a 
periodic orbit of a self-similar Sierpinski carpet billiard 0(5a); additional results 
and proofs are forthcoming in [CheNie], but a synopsis is provided in §§4.4 and 5.4. 



^The references [GaStVo, Gutl, MasTa, Sm, Tal, Ta2, Vo, Zo| provide an excellent survey of 
the various topics in the field of mathematical billiards, as well as specific results pertinent to the 
theory of rational polygonal billiards and associated flat surfaces. 

■^To our knowledge, the r-frax;tal ^ was first introduced in [AcST]. 
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Many of the results in §§4 and 5 are being announced for the first time. Specif- 
ically, §§4.3 and 5.3 contain new results on the prefractal T-fractal billiard fl{^n) 
and the T-fractal billiard ^{3^) (see [LapNie6]); §§4.4 and 5.4 contain new results 
for a prefractal Sierpinski carpet billiard fl{Sa,n) and self-similar Sierpinski carpet 
billiard Cl{Sa), where a is the single underlying scaling ratio (see [CheNie]). As 
these sections constitute announcements of new results on the respective prefractal 
and fractal billiards, we will provide in future papers [LapNie4, LapNie5, LapNie6] 
detailed statements and proofs of the results given therein. Given the nature of 
the subject of fractal billiards, we will close with a discussion of open problems and 
possible directions for future work, some of which are to appear in [CheNie] and 
[LapNie4, LapNieS, LapNie6]. 



2. Rational billiards 

In this section, we will survey the dynamical properties of a billiard ball as it 
traverses a region in the plane bounded by a closed and connected polygon. In the 
latter part of this article, we will remove the stipulation that the boundary be a 
polygon and focus on billiard tables having boundaries that are fractal or containing 
subsets that are fractal (while still being simple, closed and connected curves in the 
plane). 

Under ideal conditions, we know that a point mass making a perfectly elastic 
collision with a surface (or curve) will reflect at an angle which is equal to the 
angle of incidence, this being referred to as the law of reflection. 

Consider a compact region fl{D) in the plane with simple, closed and connected 
boundary D. Then, D,{D) is called a planar billiard when D is smooth enough to 
allow the law of reflection to hold, off of a set of measure zero (where the measure 
is taken to be the Hausdorff measure or the arc length measure on D). Though 
the law of reflection implicitly states that the angles of incidence and reflection be 
determined with respect to the normal to the line tangent at the basepoint, we 
adhere to the equivalent convention in the field of mathematical billiards that the 
vector describing the position and velocity of the billiard ball (which amounts to the 
position and angle, since we are assuming unit speed) be reflected in the tangent to 
the point of incidence. That is, employing such a law in order to determine the path 
on which the billiard ball departs after impact essentially amounts to identifying 
certain vectors. Such an equivalence relation is denoted by ~ and, in the context 
of a polygonal billiard, is discussed below in more detail. 

For the remainder of the article, unless otherwise indicated, when D is a simple, 
closed, connected and piecewise smooth curve so as to allow the law of reflection to 
hold, we assume D is a closed and connected polygon. In such a case, we will refer 
to n{D) clS cL polygonal billiard. 

We can rigorously reformulate the law of reflection, in the context of a polygonal 
billiard, as follows: the vector describing the direction of motion is the reflection — 
through the tangent at the point of collision — of the translation of the vector pre- 
viously describing the direction of motion. One may express the law of reflection in 
terms of equivalence classes of vectors by identifying these two vectors to form an 
equivalence class of vectors in the unit tangent bundle corresponding to the billiard 
table ^{D); see Figure 1. (See [Sm] for a detailed discussion of this equivalence 
relation on the unit tangent bundle fl{D) x S^.) 
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Denote by the unit circle, which we let represent all the possible directions 
(or angles) in which a billiard ball may initially move. To clearly understand how 
one forms equivalence classes from elements of fl{D) x S^, we let {x,6), {y,j) G 
0(D) X 5*^ and say that (a;, 9) {y, 7) if and only if a; = y and one of the following 
is true: 

(1) X — y is not a vertex of the boundary D and ^^ = 7; 

(2) X = y is not a vertex of the boundary D, but x = y is a point on a segment Sj 
of the polygon D and 6 = ri{-f), where denotes reflection in the segment 

'•) 

(3) If X = y is a vertex of D, then we identify {x, 9) with {y,g{'y)) for every g 
in the group generated by reflections in the two adjacent sides having x (or 
y) as a common vertex. 

For now, we shall denote by [x,9] the equivalence class of {x,9), relative to the 
equivalence relation ^. 

The collection of vertices of ^1{D) forms a set of zero measure (when we take 
our measure to be the Hausdorff measure or simply, the arc-length measure on D), 
since there are finitely many vertices. 

The phase space for the billiard dynamics is given by the quotient space (0(D) x 
S^)/ '~. In practice, one restricts his or her attention to the space {D x S^)/ '~. 
The billiard flow on {D x S^) / ^ is determined from the continuous flow on 0(_D) x 
S^)/ ~ as follows. Let a;° be an initial basepoint, 0° be an initial direction and 
ipt{x^, 9^) be a flow line corresponding to these initial conditions in the phase space 
(0(D) X S'l)/ The values for which '^t^{x",9") e {D x S'^)/ ^ constitute the 
return times (i.e., times at which ipt{x'^, 0°) returns to the section, or intersects it in 
a non-tangential way). Then, the discrete map f*^{x^,9^) constitutes the section 
map. In terms of the configuration space, p'(x^,9^) constitutes the point and 
angle of incidence in the boundary D. Since 0(D) is the billiard table and we are 
interested in determining the collision points, it is only fitting that such a map be 
called the billiard map. More succinctly, we denote by and, in general, such 
a map is called the Poincare map and the section is called the Poincare section. 
Furthermore, the obvious benefit of having a visual representation of f^{x^,9^) 
in the configuration space is exactly why one restricts his or her attention to the 
section (D x S^)/ ~. Specifically, all one really cares about in the end, from the 
perspective of studying a planar billiard, are the collision points, which are clearly 
determined by the billiard map. 

In order to understand how one determines the next collision point and direction 
of travel, we must further discuss the billiard map fn- As previously discussed, 
fn ■ {D X S^)/ ~— >■ {D X S^)/ ^, where the equivalence relation ^ is the one 
introduced above. More precisely, if 0° is an inward pointing vector at a basepoint 
x^, then {x^,9^) is the representative element of the equivalence class [(x°,0°)]. 
The billiard map then acts on (D x 5*^)/ ^ by mapping [(x'^,6''^)] to [(x'^+^, 6''^+^)], 
where x*^ and x'^^^ are collinear in the direction determined by 9'' and where 9'^'^^ 
is the reflection of angle 9'' through the tangent at x*^"*"^. In general, we have 
f^[{x°,9°)] = [(x^6l'=)], for every A; > 0. 

Remark 2.1. In the sequel, we will simply refer to an element [(.x'^, 9'')] G (0(D) x 
S^) ~ by (x'^,^'^), since the vector corresponding to 9'' is inward pointing at the 
basepoint x*^. So as not to introduce unnecessary notation, when we discuss the 
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Figure 1. A billiard ball traverses the interior of a billiard and 
collides with the boundary. The velocity vector is pointed outward 
at the point of collision. The resulting direction of flow is found by 
either reflecting the vector through the tangent or by reflecting the 
incidence vector through the normal and reversing the direction of 
the vector. We use the former method in this paper. 

billiard map corresponding to the nth prefractal billiard ri(i^„) approximating 
a fractal billiard we will simply write /f„ as /„. When discussing the discrete 

billiard flow on (f7(F„) x S*!)/ ~, the fcth point in an orbit (x*^, 6''=) G (0(F„) x 5^)/ ~ 
will instead be denoted by (a;^",6'^"), in order to keep track of the space such a 
point belongs to (namely, with our present convention, 0(F„) x S"^). Specifically, 
kn refers to the number of iterates of the billiard map /„ necessary to produce the 
pair (a:^",0^"). An initial condition of an orbit of ri(F„) will always be referred to 

as(4,0- 

An orbit with a finite period is called a closed orbit. If there exists to > 1 such 
that f^{x^,9^) = then the resulting orbit is called periodic; the smallest 

positive integer m such that f'^{x^ ,9^) = (x", 6'") is called the period of the periodic 
orbit. In the event that a basepoint of f'']^{x^ ,9^) is a corner of ^{D) (that is, 
a vertex of the polygonal boundary D), then the resulting closed orbit is said to 
be singular. In addition, if there exists a positive integer k such that the basepoint 
x~^ of /^'"'(a;°, 0°) is a corner of Vt{D) (here, denotes the fcth inverse iterate of 
/d), then the path traced out by the billiard ball connecting x^ and x~^ is called 
a saddle connection. 

We say that an orbit ff{x^,9^) is dense in a rational billiard table Q,{D) if the 
path traversed by the billiard ball in Q,{D) is dense in Q,{D). That is, the closure of 
the set of points comprising the path traversed by the billiard ball is exactly Q,{D). 
Likewise, the points of incidence (i.e., the footprint) of a dense orbit will be dense 
in the boundary D. 

Definition 2.2 (Footprint of an orbit). Let ff^ix^ ,9^) be an orbit of a polygonal 
billiard VL{D) with an initial condition (a;°,6'°) G D x 5'^ Then 



the trace of the orbit on the boundary _D, is called the footprint of the orbit 
d{x^ ,9^). When we are only interested in a prefractal billiard f2(i^„), we de- 
note the footprint of an orbit by Fn{x^,9^). 

For the remainder of the article, when discussing polygonal billiards, we will focus 
our attention on what are called rational polygonal billiards, or, more succinctly, 
rational billiards. 



(1) 
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Definition 2.3 (Rational polygon and rational billiard). If D is a nontrivial con- 
nected polygon such that for each interior angle 9j of D there are relatively prime 
integers Pj > 1 and qj > 1 such that 6j = ^w, then we call D a rational polygon 
and fl{D) a rational billiard. 

2.1. Flat surfaces and properties of the flow. In this subsection, wc will dis- 
cuss what constitutes a flat surface and how to construct a flat surface from a 
rational billiard. Then, in §2.2, we will see how to relate the continuous billiard 
flow on {Cl{D) X S^)/ with the geodesic flow on the associated flat surface. 

Definition 2.4 (Flat structure and flat surface). Let M be a compact, connected, 

orientable surface. A flat structure on M is an atlas co, consisting of charts of the 
form {Ua, Va)a6^) where Ua is a domain (i.e., a connected open set) in M and (fa 
is a homeomorphism from Ua to a domain in M?, such that the following conditions 
hold: 

(1) The collection {Ua}aes^ covers the whole surface M except for finitely 
many points z\, Z2, Zk, called singular points; 

(2) all coordinate changing functions are translations in R^; 

(3) the atlas w is maximal with respect to properties (1) and (2); 

(4) for each singular point zj , there is a positive integer mj , a punctured neigh- 
borhood Uj of Zj not containing other singular points, and a map ipj from 
this neighborhood to a punctured neighborhood Vj of a point in that is 
a shift in the local coordinates from w, and is such that each point in Vj 
has exactly mj preimages under tpj. 

We say that a connected, compact surface equipped with a flat structure is a flat 

surface. 

Remark 2.5. Note that in the literature on billiards and dynamical systems, the 
terminology and definitions pertaining to this topic are not completely uniform; see, 
for example, [GaStVo, Gutl, GutJul, GutJu2, HuSc, Mas, MasTa, Vel, Ve2, Vo, 
Zo]. We have adopted the above definition for clarity and the reader's convenience. 

We now discuss how to construct a flat surface from a rational billiard. Consider 
a rational polygonal billiard i}{D) with k sides and interior angles ^tt at each 
vertex Zj, for 1 < j < k, where the positive integers pj and qj arc relatively prime. 
The linear portions of the planar symmetries generated by reflection in the sides of 
the polygonal billiard Q{D) generate a dihedral group !^n, where N := \cm{qj}'-^^ 
(the least common multiple of the g^'s). Next, we consider f2(Z)) x H^n (equipped 
with the product topology). We want to glue 'sides' of 0(D) x together and 
construct a natural atlas on the resulting surface M so that M becomes a translate 
surface. 

As a result of the identification, the points of M that correspond to the vertices 
of J1(D) constitute (removable or nonremovable) conic singularities of the surface. 
Heuristically, ^1{D) x i^^r can be represented as {rjO(_D)}|^j^, in which case it is 
easy to see what sides are made equivalent under the action of ~. That is, ~ 
identifies opposite and parallel sides in a manner which preserves the orientation. 
See Example 2.6 and Figure 2 for an example of a flat surface constructed from the 
equilateral triangle billiard f2(A). 



FRACTAL BILLIARDS 



7 




Figure 2. The equilateral triangle billiard il{A) can be acted on 
by a particular group of symmetries to produce a flat surface that 
is topologically equivalent to the flat torus. In this flgure, we see 
that opposite and parallel sides are identified in such a way that 
the orientation is preserved. This allows us to examine the geodesic 
flow on the surface. We will see in §2.2 that the geodesic flow on 
the flat surface is dynamically equivalent to the continuous billiard 
flow. 




Figure 3. Partially unfolding an orbit of the square billiard fl{Q). 



Example 2.6. Consider the equilateral triangle A. The corresponding billiard is 
denoted by il(A). The interior angles are {■!, ^, ■!}. Hence, the group acting on 
r2(A) to produce the flat surface is the dihedral group ^3. The resulting flat surface 
is topologically equivalent to the flat torus. We will make use of this fact in the 
sequel. 

2.2. Unfolding a billiard orbit and equivalence of flows. Consider a rational 
polygonal billiard ^{D) and an orbit ff{x^^9^). Reflecting the billiard Vt{D) and 
the orbit in the side of the billiard table containing the basepoint of the orbit (or 
an element of the footprint of the orbit) partially unfolds the orbit ff[x^,9^)\ see 
Figure 3 for the case of the square billiard. Continuing this process until the orbit 
is a straight line produces as many copies of the billiard table as there are elements 
of the footprint; see Figure 4. That is, if the period of an orbit ff{x^ , 0^) is some 
positive integer p, then the number of copies of the billiard table in the unfolding 
is also p. We refer to such a straight line as the unfolding of the billiard orbit. 

Given that a rational billiard Vl{D) can be acted on by a dihedral group 3!^ to 
produce a flat surface in a way that is similar to unfolding the billiard table, we 
can quickly see how the billiard flow is dynamically equivalent to the geodesic flow; 
see Figure 5 and the corresponding caption. 

One may modify the notion of "reflecting" so as to determine orbits of billiard 
tables tiled by a rational polygon D. As an example, we consider the unit-square 
billiard table. An appropriately scaled copy of the unit-square billiard table can 
be tiled by the unit-square billiard table by making successive reflections in the 
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Figure 4. Unfolding an orbit of the square billiard 
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Figure 5 . Rearranging the unfolded copies of the unit square from 
Figure 4 and correctly identifying sides so as to recover the flat 
torus, we see that the unfolded orbit corresponds to a closed geo- 
desic of the flat surface. 



sides of the unit square. One may then unfold an orbit of the unit-square billiard 
table into a larger square billiard table. When the unfolded orbit of the original 
unit-square billiard intersects the boundary of the appropriately scaled (and larger) 
square, then one continues unfolding the billiard orbit in the direction determined 
by the law of reflection.^ We will refer to such an unfolding as a reflected-unfolding. 

We may continue this process in order to form an orbit of a larger scaled square 
billiard table. Suppose that an orbit has period p. The footprint of the 

orbit is then ^3(2:", 6*") = 0'^)}fjQ . If s is a positive integer (i.e., s G N), 

then the footprint ^|(x°,6'°) {/|j(a;°, 6'°)} -i^o"^' of an orbit constitutes the 
footprint of an orbit that traverses the same path s-many times. For sufficiently 
large s € N, an orbit that traverses the same path as an orbit (?{x^ , 9^) s-many times 
can be reflected-unfolded in an appropriately scaled square billiard table to form an 
orbit of the larger billiard table; see Figure 6. Such a tool is useful in understanding 
the relationship between the billiard flow on a rational polygonal billiard n{D) and 
a billiard table tiled by D, and will be particularly useful in understanding the 
nature of particular orbits of a self-similar Sierpinski carpet billiard in §5.4. 

As one may expect, if D is a rational polygon tiling a billiard table il{R), then an 
orbit of r2(i?) may be folded up to form an orbit of n(D). This is done by making 
successive reflections in D, the result being an orbit of r2(D); see Figure 7 for the 
case of a square billiard table. 



'That is, assuming the unfolded orbit is long enough to reach a side of the larger square. 
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Figure 6. Unfolding the orbit of the unit-square bihiard in a 
(larger) scaled copy of the unit-square billiard. This constitutes an 
example of a reflected-unfolding. The edges of the original unit- 
square billiard table and the segments comprising the orbit have 
been thickened to provide the reader with a frame of reference. 




Figure 7. Illustrated in this figure is the process of folding up 
an orbit of a square billiard table, as discussed at the end of §2.2. 
In the first image, we see an orbit of unit-square billiard table. 
Partitioning the unit square into nine equally sized squares, we see 
that we can fold up the orbit by making successive reflections in 
the sides of the squares comprising the partition. Using sufficiently 
many reflections results in an orbit of one of the squares of the 
partition. 



3. The fractals of interest 

We are primarily interested in fractals with boundaries either partially or com- 
pletely comprised of self-similar sets and fractals that are self-similar. So as to make 
the material discussed in §4-5 more accessible, we provide a few basic definitions 
from the subject of fractal geometry. 

Definition 3.1. Let {X,d) be a metric space and (p : X X. 

(i) (Contraction) . If there exists < c < 1 such that 

d{(j){x),(j){y)) < cd{x,y) 

for every x,y & X, then is called a contraction (or contraction mapping). 

(ii) (Similarity contraction). If there exists < c < 1 such that 

d{(f>{x),(t){y)) = cd{x,y), 

for every x,y X, then (j) is called a similarity contraction. This unique 
value c G (0, 1) is called the scaling ratio of (f>. 

Definition 3.2. Let {X,d) be a complete metric space. 
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(i) (Iterated function system and attractor). Let {(/>i}f=i be a family of con- 
tractions defined on X. Then {^i}'-^^ is called an iterated function system 
(IFS). ^ 

An iterated function system is so named because the map $ : K — > 
K, given by $(•) := lJi=i ^'^^ defined on the space K of nonempty 
compact subsets of X, can be composed with itself. Indeed, for each m e N, 
we have 

fc fc 

(2) U - U ^n°---°<t>iA-)- 

H = l «Tn=l 

Furthermore, there exists a unique nonempty compact set F C X (i.e., 
F gK), called the attractor of the IFS, such that 

fc 

(3) F = ^F):=\JUF)- 

i=l 

(ii) (Self-similar system and self-similar set). In the special case where each 
(f)i is a contraction similarity, for i = 1, k, then the IFS {<j)i}^^l is said to 
be a self-similar system and its attractor F is called a self-similar set (or 
a self-similar subset of X). 

If X is complete, then so is K (equipped with the Hausdorff metric) and hence, 
since it can be shown that $ : K — > K is a contraction, it follows from the con- 
traction mapping theorem that $ has a unique fixed point (thereby justifying the 
definition of the attractor F above) and that for any £^ € K, $™(i?) F, as 
m — )■ oo (where, as in Equation (2), $™ is the mth iterate of $). (See [Hut].) 

We state the next property in the special case that will be of interest to us, 
namely, that of an IFS in a Euclidean space. 

Theorem 3.3 ([Hut]; see also [Fa, Thm. 9.1]). Consider an iterated function 
system given by contractions {<j>i}^^i. each defined on a, compact set D C M", such 
that 4>i{D) ^ D for each i < k, and with attractor F. Then F C D and in fact, 

oo 

(4) ^ = n '^"'(E) 

m=0 

for every set E G K such that (j)i{E) C E for all i < k. Here, the transformation 
$ : K ^ K is given as in part [i) of Definition 3.2. 

Not every fractal is self-similar or embedded in Euclidean space. However, such 
sets represent an important collection of examples of fractal sets. In the next 
subsection, we will discuss the fractal subsets (self-similar or not) of M or of of 
direct interest to us in this paper. 

3.1. Cantor sets. A Cantor set is a set with very rich and counter-intuitive prop- 
erties; topologically, it is a compact and totally disconnected (i.e., perfect) space. 
In order to illustrate some of the properties that make a Cantor set so interesting, 
we refer to the canonical example of a Cantor set: the ternary Cantor set. We focus 
on three methods for constructing the ternary Cantor set: 1) by tremas, 2) as the 
unique fixed point attractor of an iterated function system, and 3) in terms of an 
alphabet. 
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Figure 8. The ternary Cantor set. 



Before we discuss the ternary Cantor set, we comment on the fact that such 
a set was first discovered by Henry J. S. Smith in 1875. Later, in 1881, Vito 
Voltcrra independently rediscovered the ternary Cantor set. Smith's and Volterra's 
records being obscured over the years for one reason or another, it was the German 
mathematician Georg Cantor whom, in 1883, history credits with the discovery of a 
bounded, totally disconnected, perfect and uncountable set with measure zero that 
is now commonly referred to as "the Cantor set". 

We now proceed to construct the ternary Cantor set, hereafter denoted by , 
by the method known as construction by tremas, which is Latin for 'cuts'. Begin 
with the unit interval / and remove the middle open third (5, §) from /, leaving 
the two closed intervals [0, |] and [|, 1]. Next, remove the middle open ninth from 
each closed subinterval. What remains are the closed intervals [0, [§,5], [|, |], 
[|, 1]. Continuing this process ad infinitum, we construct the ternary Cantor set; 
see Figure 8. 

One may also construct 'rf set by utilizing an appropriately defined iterated 
function system. Consider the following contraction maps on the real line M: 

1 12 
(5) = 3 a;, =3^+3- 

Then, lim„_,.oo = where $ is given as in part (i) of Definition 3.2. More- 

over, since {4)i}f^i is a family of similarity contractions and = $(^), we have 
that is a self-similar set. 

A third — and equivalent — construction of the ternary Cantor set can be given 
in terms of the symbols I, c, and r. Recall that elements of R can be expressed in 
terms of a base-3 number system. We focus our attention on elements of the unit 
interval I. So-called ternary numbers* in I have two equivalent expansions: one 
that is finite and one that is infinite. For example, i can be written in base-3 as 0.1 
or, equivalently, as 0.02 (where the over bar indicates that the digit 2 is repeated 
infinitely often). 

We next discuss a similar addressing system that has the benefit of preventing 
ternary numbers from having a finite representation. The characters I, c and r are 
to remind the reader of choosing left, center and right. We identify an element of 
the unit interval / by an infinite address that indicates where in / the element is 
located. Motivated by the construction of 'rf by tremas, one can identify any point of 
I by an infinite address consisting of the characters I, c and r. While elements of ^ 
can be represented by infinite addresses consisting of c's, we make the stipulation 
that no element of "if will be represented by an infinite address containing c's.^ 



An element x' G 7 is a ternary number if a; = < p < 3^, p, j/ G N. 
^In other words, we do not allow an element of to be approximated by a sequence {zi}°. 
of elements of where is the complement of the ternary Cantor set in I. 
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Moreover, this method of representing elements of / (or 'rf) provides every element 
with an infinite representation and never a finite representation. 

Example 3.4. The values j, ^ and ^ have the ternary representations Ir, If and 

c, respectively.® While i has a finite ternary expansion given by 0.2 = 0.01, it docs 
not have a finite ternary representation. It should be noted that elements like ^ and 
^ will play an important role in our analysis of the Koch snowflake fractal billiard. 

The occurrence of infinitely many c's or infinitely many Vs and r's is critical to 
developing some of the theory regarding the Koch snowflake fractal billiard. 

So that some of the results concerning the Koch snowflake fractal billiard can be 
more succinctly expressed, we define a notation used for describing a value's type 
of ternary representation. 

Notation 3.5 (The type of ternary representation). The type of ternary repre- 
sentation can be defined as follows. If a: e /, then the first coordinate of [•,•] 
describes the characters that occur infinitely often and the second coordinate of 
[•, •] describes the characters that occur finitely often. If we want to discuss many 
different types of ternary representations, then we use 'or'. That is, the notation 
[•, •] V [•, •] V ... V [•, •] is to be read as [•, •] or [•, •] or ... or [•, •]. If the collection of 
characters occurring finitely often is empty, then we denote the corresponding type 
of ternary representation by [-,0]. 

Example 3.6. The value \ has a ternary representation of c. Hence, \ has a type 
of ternary representation given by [c, 0] . Moreover, the value ^ has a ternary rep- 
resentation given by crl, which means that ^ has a type of ternary representation 
given by [Zr, c]. 

We note that "the" type of representation of a point .t G J is not unique, in 
general. For instance, the value ^ has a ternary representation of [r,l] or [l,c\. 

A thorough understanding of the ternary Cantor set is not only important for 
understanding many of the results on the Koch snowflake prefractal and fractal 
billiard. In general. Cantor sets will be ubiquitous and instrumental in our anal- 
ysis of other fractal billiard tables. In each example of a fractal billiard, we will 
clearly indicate where and how a particular Cantor set is important in analyzing a 
particular fractal billiard table. 

3.2. The Koch curve and Koch snowflake. The Koch curve KC is constructed 
as shown in Figure 9 and is the unique fixed point attractor of the following iterated 
function system on the Euclidean plane (here, i = \/— T): 

(6) <^i(x) = ^x, ^2(x) = ie^?x+(i,0), 

^3(x) = ie-^f X + (| ^), </.4(x) = ix + (| 0). 

Since each contraction map in the iterated function system is a similarity transfor- 
mation (i.e., {<t>j}j=i is a self-similar system) and KC = ^{KC), we have that KC 
is a self-similar set; see part (ii) of Definition 3.2. There are additional properties 



Equivalently, g has a representation given by cl. Although, we will not consider this as a 
representation for g on account of g G 
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Figure 9. The construction of the Koch curve KC. Here, the self- 
similar set KC is viewed as a limit of the prefractal approximations 
{$'"(i?)}^^0) where R is the initial triangle and the map i> is 
defined in terms of the IFS given by Equation (6), as in Definition 
3.2. (See Theorem 3.3 and the text preceding it.) 



Figure 10. One typically sees this construction of the Koch curve 
KC when learning about fractal sets. Beginning with the unit 
interval I, one removes the middle third and replaces it with the 
two other sides of an equilateral triangle, as shown. One then 
repeats this process infinitely often for every remaining interval; 
the resulting limiting set is KC. Such a sequence {i^C„}5^g of 
approximations converges to KC because it is a subsequence of 
the convergent sequence of prefractal approximations {$'"(i?)}^^Q 
shown in Figure 9. (Here, we are using the notion of convergence 
in the sense of the Hausdorff metric.) 



of the Koch curve that are reminiscent of the Cantor set, which is more than just 
a coincidence, this being discussed in more detail below. 

If we allow the iterated function system to act on the triangle R = {{x,y)\0 < 
X < ^,0 < y < ^x} U {{x,y)\^ < x < 1,0 < y < -^x + ^}, as shown in 
Figure 9, sequential iterates of the iterated function system very quickly produce 
a prefractal that is visually indiscernible from the true limiting set. But there is a 
more common construction that allows us to visualize the curve KC more readily, 
this being depicted in Figure 10. The technical caveat which we are brushing under 
the carpet is that each polygonal approximation shown in Figure 10 does not contain 
the Koch curve KC, while each approximation in the sequence shown in Figure 9 
does contain KC."^ 

Notation 3.7. For each integer n > 0, we denote by iS7„ and KSn the nth 
(inner) polygonal approximations of the Koch curve KC and Koch snowflake KS, 
respectively. 



Recall from Definition 3.2 and Theorem 3.3 that for a set -F to be the unique fixed point 
attractor of an IFS, each F„ must be such that F C $(Fn), so that $(Fn) = *"+^(-Fb). 



A 




KCo 



KC 
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Figure 11. The Koch snowflake is comprised of three Koch 
curves. We have encapsulated each Koch curve to highlight how 
KS is the union of three abutting copies of KC . 



Intuitively, one expects the Koch curve to have finite length, since it is the limit 
of a sequence of polygonal approximations. On the contrary, the Koch curve KC has 
infinite length, which can be seen by the following calculation given in terms of the 
nth prefractal KC'm where KCn is one of the polygonal approximations indicated 
in Figure 10: 

(7) length of KCn = 

Then, lim„^oo (|)" = oo. 

The Koch snowflake KS is a fractal comprised of three abutting copies of the self- 
similar Koch curve; see Figure 11. As a closed (simple) curve, the Koch snowflake 
KS bounds a region of the plane; furthermore, the area of this region can be cal- 
culated as follows: 

n 

(8) area bounded by KSn = 1 + 

i=0 

Then, as n increases, the right-hand side of (8) tends to a finite value. The area 

bounded by the Koch snowfiake is thus given by lim„_i.oo 1 + S"=o (f )* ~ ^' 
suming the sides of KSq have length one. 

As we noted at the end of §3.1, Cantor sets are ubiquitous in the context of 
self-similarity. In the case of the Koch snowflake, KS n KSn is the union of 3 • 4" 
self-similar ternary Cantor sets, each spanning a distance of Such a fact will be 
important in determining certain sequences of what we will call compatible orbits 
(see Definitions 4.2-4.5) and certain families of well-defined orbits of il{KS). 

3.3. The T-fractal. The T-fractal ^, first introduced in [AcST], is not a self- 
similar set. However, contains a Cantor set and is constructed in a way that 
is reminiscent of an IFS acting on a compact set to produce a self-similar set.® 
As shown in Figure 12, one constructs the T-fractal by appending scaled copies 



'^Recall from Definition 3.2 and Tlieorem 3.3 tiiat each prefractal approximation Fn must 
contain tlie unique fixed point attractor F. 
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Figure 12. The construction of the T-fractal . 



Figure 13. can be tiled by the unit square Q. 



of the initial T shape to each successive approximation. Specifically, S/'n+i is 
constructed from .%i by appropriately appending 2"+-'^ copies of 2stt=% to S/'n^ 

The overall height of can be calculated and the total area bounded by can 
be shown to be finite, as shown in the following calculations (we assume here that 
the base of =3^ is two units in length): 

(9) heightof^„ = 3+^ + ^ + ... + |^=3f^i^. 

i=0 

Then, lim„_,.(x) 3 X^^g ^ ~ ^' which is the height of . Furthermore, the area 

bounded by S^n is calculated as follows. There are eight squares, each with side- 
length one, comprising 3^; see Figure 13. Hence, the area of 3'q is eight square- 
units. Therefore, 

(10) area bounded by = 8 2 • + 2" • — = 8 V — . 

i=0 

Then, lim„^oo 8 ^"^q ^ = 16, which is the total area bounded by ,5^. 

There is a natural fractal subset of but no point of is in this fractal subset, 
which is unlike what we have seen in the case of the Koch snowflake fractal KS. In 

fact, the fractal subset in question is given by {(x, 6)|x G R} n 5^. We note that 
this fractal subset is not self-similar, though it is a (topological) Cantor set. 

3.4. Self-similar Sierpinski carpets. A Sierpinski carpet can be constructed 
by systematically removing particular open subsquares from the unit square Q = 
{{x,y)\0 < X < l,0<y< 1}. Depending on how one chooses the sizes of the open 
subsquares to be removed, one can either construct a self-similar Sierpinski carpet 
or a non-self-similar Sierpinski carpet, these being defined below. This method 
of construction is called construction by tremas and is described in the caption 



9See [AcST, §2.1] and [LapNie6] for a more precise description of the definition of 
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Figure 14. The 1/3-Sierpinski carpet is a self-similar carpet con- 
structed in one of two ways: 1) by tremas and 2) an iterated func- 
tion system (in fact, a self-similar system). We describe here the 
construction of the 1/3-Sierpinski carpet by tremas, the latter be- 
ing discussed in the main text. Beginning with the unit square, 
one then removes the middle open square with side-length ^ . From 
each remaining subsquare of side-length i, one then removes the 
middle open square of side-length ^. One continues this procedure 
of removing subsquares of remaining squares until there is no area 
left. As one would expect, each step of the construction process 
can be emulated by applying the correct iterated function system, 
which is given in Equation (11). 



of Figure 14, using the standard "1/3-Sierpinski carpet" as an example. Such a 
construction process should be very familiar, since "removing middle thirds" is 
exactly what we did to construct the ternary Cantor set in §3.1. 

As referred to in the caption of Figure 14, one may also construct the 1/3- 
Sierpinski carpet by applying an appropriately defined iterated function system to 
the unit square Q. Consider the following iterated function system, which is a 
self-similar system. 

(11) -^iW = ^x, <^2(x) = + (0, \ 



3 ' ' 3 V ' 3 

(/.3(x) = ( 0,^ ) , 04(x) = ix+(io 



3 V 3y ' ' 3 V3 

</,,(x) = -x+(^-,-j, 0s(x) = -X+(^-,- 

Then, denoting the 1/3-Sierpinski carpet by 6*3, we have that lim„^oo $"((3) = 6*3. 
Since each contraction in the iterated function system is a similarity contraction 
and 53 = ^{Ss), it follows that ^3 is a self-similar set. 

We discuss here the relevant results and material from [Du-CaTy]. For our 
purposes, the first level approximation of a Sierpinski carpet Sg, will always be the 
unit square Q and denoted by Sq. Since every (self-similar and non-self-similar) 
Sierpinski carpet has the same zeroth level approximation and zero is never a scaling 
ratio, such notation will never cause any confusion. 

What follows is a general description on how to construct a Sierpinski carpet 
by removing appropriately sized middle, open squares. Consider the unit square 
Q = Sq. Let ao = 2fco + 1 for some fco S N. Partition Sq into ao squares of 
side-length a^^. Next, remove the middle, open subsquare. Let a\ = 2ki + 1 for 
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some fci G N. Each subsquare may then be partitioned into many squares with 
side-length (ao • ai)~^- We then remove each middle, open subsquare of side-length 
(ao • ai)~^; see Figure 14. Continuing this process, let a„_i = 2A:„_i + 1 where 
kn-i € N and let a„ = 2fc„ + 1 for some fc„ e A''. Then we partition a subsquare of 
side-length (ao • ai • • • a„_i)~^ into many squares. We then remove the middle, 
open square from each subsquare in the partition. Continuing in this manner ad 
infinitum, one constructs a Sierpinski carpet denoted by ^a, where a = {a~^}^g. 

Definition 3.8 (A self-similar Sierpinski carpet). If a = {a^^}iZoy with a, = 2ki+l 
and fcj e N, is a periodic sequence of rational values, then the Sierpinski carpet 5a 
is called a self-similar Sierpinski carpet. 

We have described the construction of a self-similar Sierpinski carpet Sa. in terms 
of the removal of particular open squares. As the name would suggest, there exists 
an appropriately defined iterated function system such that S'a = $(<S'a). 

Viewing Sa. as the unique fixed point attractor of an appropriately defined iterated 
function system will be useful in proving some of the statements in the subsequent 
sections. More precisely, Sa is viewed as the self-similar set associated with a self- 
similar system, as in part (ii) of Definition 3.2. 

While we do not discuss any non-self-similar Sierpinski carpet billiards in this 
paper, we provide the definition for completeness. 

Definition 3.9 (A non-self-similar Sierpinski carpet). If a = {a~^}^Q, with a, = 
2ki + 1 and ki £ N, is an aperiodic sequence of rational values, then the Sierpinski 
carpet 5a is called a non-self-similar Sierpinski carpet. 

Definition 3.10 (A cell of 5a.n)- Let ao = 2fco -{- I, ko E N. Consider a partition 
of the unit square Q — So into many squares of side- length a^^. A subsquare of 
the partition is called a cell of So and is denoted by Co,ao ■ Let 5a be a Sierpinski 
carpet. Consider a partition of the prefractal approximation 5a,n into subsquares 
with side-length (ao • ai • • • a„)~^. A subsquare of the partition of 5a,n is called a 
cell of Sa,n and is denoted by Cn,aoai---an and has side-length (ao • ai • • • a„)~^. 

Definition 3.11 (Peripheral square). In accordance with the convention estab- 
lished in [Du-CaTy], the boundary of an open square removed in the construction 
of 5a is called a peripheral square of 5a. Furthermore, by convention, the unit 
square Q = 5o is not a peripheral square. 

Definition 3.12 (Nontrivial line segment of 5a). A nontrivial line segment o/5a 
is a (straight-line) segment of the plane contained in 5a that has nonzero length. 

Unless otherwise indicated, in what follows, we assume that 5a is a self-similar 
Sierpinski carpet with a single scaling ratio a; that is, a = {a~^}^Q, where a = 
2k +1 for some fixed fc G N. In addition, when a = {a~^}°Zo^ '^a is denoted by Sa- 

We next state the following theorem, due to Durane-Cartagena and Tyson in 
[Du-CaTy] and which will be very useful to us in this context (see §4.4 and §5.4). 

Theorem 3.13 ([Du-CaTy, Thm. 4.1]). Let Sa be a self-similar Sierpinski carpet. 
Then the set of slopes Slope{Sa) of nontrivial line segments of Sa is the union of 
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the following two sets: 

(12) A = : p + q < a, 0<p<q<a-l, e NU{0}, p + q is odd 

(13) B = 1^ : p + g < a - 1, < p < g < a - 2, p,qeN, p,q are odd 

Moreover, if a G A. then each nontrivial line segment in Sa with slope a touches 
vertices of peripheral squares, while if a £ B, then each nontrivial line segment in 
Sa with slope a is disjoint from all peripheral squares. 

Notation 3.14. Let a,b be odd positive integers such that 3 < 6 < o and let 
Slopc(S'a) and Slope(S'f,) be the set of slopes of nontrivial line segments of Sa and 
Sb, respectively. We denote by Aa (resp., Ab) the subset A C Slope(S'a) (resp., 
A C Slope(S'h)) given in Equation (12) of Theorem 3.13. Similarly, we denote by 
Ba (resp., Bb) the subset B C Slope(S'o) (resp., Slope(S'6)) given in Equation (13) 
of Theorem 3.13.^" 

If Sa and Sb are self-similar Sierpinski carpets with b < a, then it is clear that 
Slope(S'6) C Slope(S'(i). Moreover, in this case, we also have that A^ C Aa and 
Bb C Ba. 



4. Prefractal (rational) billiards 

111 the previous sections, we surveyed basic facts and results from mathematical 
billiards and fractal geometry, with most of our attention being focused on the 
subject of rational billiards and sets exhibiting self-similarity. We also discussed the 
importance of examining the dynamically equivalent geodesic flow on an associated 
flat surface. In this section, we will examine examples from particular classes of 
prefractal (rational) billiards. We are interested in tables that can be tiled by a 
single polygon that can also tile the (Euclidean) plane. The main examples we 
will discuss are the Koch snowflake prefractal billiard table, the T-fractal prefractal 
billiard table and a self-similar Sierpinski carpet prefractal billiard table. Each 
example of a prefractal billiard table constitutes a rational billiard table, but is 
an element of a sequence of rational billiard tables approximating a fractal billiard 
table with radically diS'erent qualities when compared to the others. That is, the 
Koch snowflake has an everywhere nondiflferentiable boundary; the T-fractal billiard 
table is certainly a fractal billiard table, since its boundary ^ contains a fractal 
set, but the portion of the boundary that is nondifferentiable has Lebesgue measure 
zero; a Sierpinski carpet billiard table can possibly have no area, yet yield billiard 
orbits of flnite length. 

4.1. A general structure. We restrict our attention to billiard tables with fractal 
boundary F, where F can be approximated by a suitably chosen sequence of rational 
polygons {Fn}^^Q- More specifically, we are interested in a fractal billiard table 
fl{F) with the property that, for every n > 0, ri(F„) can be tiled by a single polygon 
Dn, where Z)„ = c„£>o- Here, < c„ < 1 is a suitably chosen scaling ratio and Dq 



In the case of Af, (resp., Bf,), a should of course be replaced by b in Equation (12) (resp., 
Equation (13)). 
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Figure 15. In the first image, we liave tlie orbit i^o(c, 7r/3) of 
il(A'5'o). In the second image, we see that the orbit ^o{c, f), when 
embedded in rt{KSi), is not an orbit ofil.{KSi). In the third image, 
the orbit intersects the sides of n{KSi) and appears to be related 
to (c, f ) in some way. 



is a polygon that tiles both the (Euclidean) plane as well as the rational billiard 

The focus in this subsection is on developing a general framework for discussing 
billiards on prefractal approximations. If n{Fn) and f2(i^„-|-i) are two prefractal 
billiard tables approximating a given fractal billiard table ^{F), then we want 
to have a systematic way of determining how and if two orbits ^„(a;^,^J^) and 
^«+i(a;°+i, 6*^+1) of n{Fn) and n{F„+i) are related. 

Notation 4.1. We will primarily measure angles relative to a fixed coordinate 
system, with the origin being fixed at a corner of a prefractal approximation Fq. 
However, we will sometimes measure an angle relative to a side of Fn on which a 
billiard ball lies. In such situations, we will write the angle as w{9) in order to 
indicate that the inward pointing direction is 9, measured relative to the side on 
which the vector is based. 

To motivate our general discussion, consider the orbit ^0(2^0, 7''/3) of il.{KSo), 
where Xq = c € /; see the first image in Figure 15 (and recall our earlier discussion in 
§3.1). The same orbit, viewed as a continuous curve embedded in fl(KSi), does not 
constitute an orbit of n{KSi); see the second image in Figure 15. Consider the orbit 
^i{xi, ^) shown in the third image in Figure 15. Such an orbit does intersect the 
boundary of fl(KSi) and appears to be related to &o(c, but in what way we have 
not yet explicitly said. Initially, we notice that, as a continuous curve embedded 
in il{KSi), the orbit ^o{c, -1) is a subset of ffi{xi, Being eager to establish a 
proper notion of "related", we may be inclined to declare that two orbits are related 
if one is a subset of the other, when viewed as continuous curves in the plane. 
Unfortunately, we quickly see that such a definition is highly restrictive. A more 
general observation is that and Xq are coUinear in the direction of ^, without 

any portion of KSi intersecting the segment x^Xq. We then say that {xq, ^) and 
{xi, ^) are compatible initial conditions. We state the formal definition as follows. 

Definition 4.2 (Compatible initial conditions). Without loss of generality, suppose 
that n and m are nonnegative integers such that n > m. Let (x^, 9^) e {Fn xS^)/ ^ 
and (x^, 9^) £ {Fm x S*^)/ ^ be two initial conditions of the orbits i^?'„(x" , 9^) and 

^^In the case of certain prefractal approximations, n(Fo) is exactly Dq. In general, however, 
Do does not always equal f2(Fo), but certainly tiles f2(Fo). An example of this situation is f!(,ifo); 
see §4.3. Such a billiard is tiled by the unit square, which is the associated polygon Dq. 
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^mC-Tmi^m); respectively, where we are assuniing that 0^ and 6'°„ are both inward 
pointing. If O'-^ = 6^^^ and if and a;^ he on a segment detorminod from 6^ (or 6*^^) 
that intersects F„ only at , then we say that (a;° , and (a;^^, 0^) are compatible 
initial conditions. 

Remark 4.3. When two initial conditions (a;°,6'^) and {x'}^,9^) are compatible, 
then we simply write each as {x^,6^) and (a;^,6'°). If two orbits i^rn{x'^,0^) 
and ffn{x'^,6'^) have compatible initial conditions, then we say such orbits are 

compatible. 

Depending on the nature of i^{F), not every orbit must pass through the region 
of ri(F„) corresponding to the interior of Q{Fo), let alone pass through the interior 
of n(Fm), for any m < n. Because of this, it may be the case that an initial 
condition is not compatible with (a;^,0°), for any m < n. As such, in 

Definitions 4.4 and 4.5, we consider sequences beginning at i = N, for some N >0. 

Definition 4.4 (Sequence of compatible initial conditions). Let {{Xi,9^)}?^ff be a 
sequence of initial conditions, for some integer iV > 0. We say that this sequence is 
a sequence of compatible initial conditions if for every m, > N and for every n > m,, 
we have that {x^,6^) and {x'^,6^) are compatible initial conditions. In such a 
case, we then write the sequence as {(a;°, 

Definition 4.5 (Sequence of compatible orbits). Consider a sequence of compat- 
ible initial conditions {(a;J^, ^°)}J^jv Then the corresponding sequence of orbits 
{^n{Xn, 6^)}^=N is called a sequence of compatible orbits. 

If ffrn{x^, 0^) is an orbit of fl{Fm), then ^^(a;^, 6'^) is a member of a sequence 
of compatible orbits {^„(a;" , 6'°)}5^^ for some > 0. It is clear from the definition 
of a sequence of compatible orbits that such a sequence is uniquely determined by 
the first orbit ffN{x%,6^). Since the initial condition of an orbit determines the 
orbit, we can say without any ambiguity that a sequence of compatible orbits is 
determined by an initial condition (x'^,6'^). 

Definition 4.6 (A sequence of compatible P orbits). Let P he a. property (resp., 
Pi,...,Pj a list of properties). If every orbit in a sequence of compatible orbits has 

the property P (resp., a list of properties Pi, ...,Pj), then we call such a sequence 
a sequence of compatible P (resp., Pi, ...,Pj) orbits. 

The following theorem can be deduced from Theorem 3 of Gutkin's paper [Gut2]; 
see [LapNie3, §3.2]. 

Theorem 4.7. Consider a prefractal rational billiard f2(F„). Ifi}{Fn) is tiled by a 
rational polygon D„ such that Dn tiles the Euclidean plane, then an orbit ^„(x^, 6^) 
ofQ,(Fn) is either closed or dense in r2(i^„)/^ regardless of the initial basepoint x^. 

Remark 4.8. When ri(i^„) is tiled by where Dn is a rational polygon tiling 
the plane, then ri(i^„) is more generally referred to as an almost integrable billiard, 
this being the language used in [Gut2]. 

The following is a generalization to this broader setting of Corollary 16 from 
[LapNie3]. It is established in the same manner. 



'Recall that these notions were introduced towards the beginning of §2. 
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Theorem 4.9. Let ^1{F) be a fractal billiard table approximated by a suitable se- 
quence of rational polygonal billiard tables {Q{Fn)}^^Q. If there exists a, polygon 
Dq that tiles the plane and such that for every n > there exists < c„ < 1 with 
Dn '■= CnDo tiling i}{Fn), then any sequence of compatible orbits is either entirely 
comprised of closed orbits or entirely comprised of orbits that are dense in their 
respective billiard tables. 

4.2. The prefractal Koch snowflake bilHard. The billiard n{KSn) can be tiled 
by equilateral triangles. Specifically, if A is the equilateral triangle with sides having 
unit length, then n{KSn) is tiled by ;^A, for every n > 0. Moreover, as is well 
known, A = KSq tiles the plane. Therefore, Theorems 4.7 and 4.9 hold for the 
prefractal billiard il{KSn)- 

Our goal for this subsection and §4.2.1 is to survey some of the main results 
of [LapNiel, LapNie2, LapNieS]. We will focus on pertinent examples that will 
motivate a richer discussion in §5.2. Initially, we focus on properties of orbits with 
an initial direction of ^ and f.^"^ 

3 6 

If ^q{xq, ^) is an orbit of il(_KS'o), so long as .Tq is not a corner of KSq, the orbit 
will be periodic, as expected. However, depending on the nature of the ternary 
representation of Xq, the compatible orbit ^i(a;5, f ) may be singular in Cl{KSi)}^ 

Theorem 4.10 ([LapNie3]). Let Xq G I C KSq. If Xq has a ternary representa- 
tion of type [l,cr] V [r,lc], then there exists N > such that the compatible orbit 
^n{x%i, ^) will be singular in Q.{KSn)- Moreover, for every n> N , ^„(a;^, |-) will 
be also singular in Sl{KSn)- 

Theorem 4.11 ([LapNieS]). If Xq has a ternary representation which is of the 

form [c,lr] V [lc,r] V [cr,l] V [^cr, 0] V [lr,c], then the sequence of compatible orbits 
{^n{x'^,0'^)}^^Q is a sequence of compatible periodic orbits. 

Theorem 4.12. The length and period of an orbit £^m{x%^, ^) G {^^(xj^, f )}J^o 
dictated by the ternary representation of Xq. {See [LapNie2] for the corresponding 
specific formulas.) 

Proof. See §4.4 of [LapNic2] for a proof of this result, as well as for additional 
properties of orbits with an initial direction of 7r/3. ■ 

Example 4.13 (A sequence of compatible hook orbits). Let Xq & I have a ternary 
representation given by rl. Such a representation indicates that, in each prefractal 
approximation KSn, a^o ''^^ element of an open, connected neighborhood con- 
tained in KSn- The point Xq corresponds to the value 3/4 e /. If we consider 
an orbit of Q.{KSq) with an initial direction of 7r/6, the ternary representation of 
the basepoints at which the billiard ball path forms right angles with the sides of 
n(iCS'o) is of the type [c,lr\. This is a degenerate periodic hybrid orbit, meaning 
that it doubles back on itself, and the next orbit in the sequence of compatible 
periodic hybrid orbits has the initial condition (.T5,7r/6) = (a;Q,7r/6). Since the 
ternary representation of the basepoint of /o(a;Q,7r/6) is rc and = = 7r/6, it 



Equivalently, we could focus on orbits with an initial direction of ^ and ^, since ^ is the 
rotation of through the angle , the angle being an angle that determines an axis of symmetry 
oi KSn, for 71 > 0. 

^'^Recall that the notion of type of a ternary representation was introduced in Notation 3.5 of 
§3.1. 
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follows that the basepoint of 7r/6) is a point which, for every prefractal ap- 

proximation KSn, is an element of an open, connected neighborhood contained in 
KSn- Then the basepoint of fi{x\, 7r/6) (where denotes the second iterate of the 
billiard map /i) has a ternary representation of type [c^lr]. This same pattern is 
repeated for every subsequent orbit in the sequence of compatible orbits. It follows 
that the resulting sequence of compatible orbits forms a sequence of orbits that is 
converging to a set which is well defined. That is, such a set will be some path in 
the fractal billiard table D,{KS) with finite length which is effectively determined 
by the law of reflection in each prefractal approximation of Q,{KS). 

Such orbits are introduced in [LapNieS] and referred to as hook orbits^ because 
they appear to be "hooking" into the Koch snowflake; see Figure 16. 

The hook orbits of Example 4.13 are special cases of a general class of orbits 
called hybrid orbits and also introduced, as well as studied, in [LapNieS]. 

Definition 4.14 (Hybrid orbit). Let ^„(a;° , 6^°) be an orbit of n{KS.„). If aU but 
at most two basepoints a;fj" e J-n{x^^ 6^) have ternary representations (determined 
with respect to the side s„^^ on which each point resides) of type [c, Zr] V [cZ,r] V 
[cr, I] V [Icr, 0] V [Ir, 0], then we call ,^„(a;° , 6lO) a hybrid orbit of n{KSn)- 

Definition 4.15 (A 3^ hybrid orbit). If ^„(a;° , 0°) is a hybrid orbit with property 
then we say that it is a ^ hybrid orbit. 

Proposition 4.16. // (x^ , 6*,",) is a dense orbit of fl(KSn), then ^„(x^j,0^) is a 
dense hybrid orbit. 

Applying the results in Theorem 4.9 and Proposition 4.16, we state the following 
result. 

Theorem 4.17 (A topological dichotomy for a sequence of compatible orbits, 
[LapNieS]). Let {^n(2^^, ^'^)}^7v sequence of compatible orbits. Then we have 
that the sequence of compatible orbits {i^„(a;° , 0'')}^^^ is either entirely comprised 
of closed orbits or is entirely comprised of dense hybrid orbits.^^ 

Theorem 4.18 ([LapNieS[). If ffo{x^,0^) is a periodic hybrid orbit of n{KSo) 
with no basepoints corresponding to ternary points {i.e., points having ternary rep- 
resentations of the types [l,cr] V [r, Zc]), then for every n> 0, the compatible orbit 
^n{x'!^,0n) ^s aperiodic hybrid orbit ofQ,{KSn)- 



Recall that the notions of "closed orbit" and "dense orbit" were defined in §2, just before 
Definition 2.2. 
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Figure 17. Three examples of periodic hybrid orbits. These are 
the first three elements of the sequence of compatible periodic hy- 
brid orbits described in Example 4.20. 



In order to fully understand the following result, we define what it means for a 
vector to be rational with respect to a basis {ui, U2} of R^. If 2; = mui + nu2, for 
some m,n € then we say that z is rational with respect to the basis {ui,U2}. 
Otherwise, we say that z is irrational with respect to {^1,7/2}. 

Theorem 4.19 (A sequence of compatible periodic hybrid orbits, [LapNieS]). 

Let Xq G I and consider a vector (a, b) that is rational with respect to the basis 
{1x1,1*2} := {(1, 0), (1/2, -v/S/S)}. Then, we have the following: 

(1) // a and b are both positive integers with b being odd, Xq — j^, for some 
r, s G N with s > 1, 1 < r < 4* being odd and 9^ := arctan^^^, then 
the sequence of compatible closed orbits {i^nla^n, ^")}^o ^ sequence of 
compatible periodic hybrid orbits. 

(2) If a = 1/2 , b is a positive odd integer, Xq = for some r, s € N with s > 1, 

1 < r < 2" being odd and 6*^ := arctan^^, then the sequence of compatible 
closed orbits {^n(a;^, ^*')}5^q is a sequence of compatible periodic hybrid 
orbits. 

Example 4.20 (A sequence of compatible periodic hybrid orbits). In Figure 17, 
three periodic hybrid orbits are displayed. These three orbits constitute the first 
three terms in a sequence of compatible periodic hybrid orbits. If we choose Xq — 
c £ I and 9q to be an angle such that Xq connects with the midpoint of the lower 
one-third interval on the side of il{KSo), we can see that ffo{xQ,6Q) is a periodic 
hybrid orbit. More importantly, there are elements of the footprint TqIxq, 9q) with 
ternary representations of type [Ir, c] . This observation is key for constructing what 
we call nontrivial paths of Q{KS), a topic which is discussed in more detail in §5.2. 

Theorem 4.21 (A constant sequence of compatible periodic hybrid orbits, 
[LapNieS]). If every element ofJ-o{xQ, 9q) has a ternary representation of type [Ir, c], 
then there exists N >0 such that {^Ti(ic" , 0*')}J^jy constant sequence of com- 
patible periodic hybrid orbits. 

Example 4.22 (A constant sequence of compatible periodic hybrid orbits). Con- 
sider Xq = 7/12 in the base of the equilateral triangle. Such a value has a ternary 
representation of type [Zr, c]. Consider the initial condition (xq, Then the se- 
quence of compatible orbits {€?„(a;° , ■|)}5^i is a constant sequence. This follows 
from the fact that the ternary representation of Xi is rl. Moreover, the represen- 
tation of every basepoint of ^„(x^, f) is Ir. In Figure 18, we show the first three 
orbits in this (eventually) constant sequence of compatible periodic hybrid orbits. 
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Figure 18. An eventually constant sequence of compatible peri- 
odic hybrid orbits. We see that the initial basepoint Xq = 7/12 lies 
on the middle third of the unit interval. The basepoint xl of the 
compatible initial condition ^) has a ternary representation of 
type [/r,0]. 




Figure 19. The flat surfaces S{KSi), S{KS2) and SiKS^) associ- 
ated with the Koch snowflake prefractal approximations KSi , KS2 
and KS3, respectively. 



As of now, the only examples of constant sequences of compatible nondegenerate 
periodic hybrid orbits we can provide are those for which the initial direction is ^. 
Of course, one can construct a constant sequence of compatible periodic hybrid or- 
bits, each with an initial direction of 7r/6 or n/2, but such orbits will be degenerate, 
as discussed in Example 4.13. 

4.2.1. The corresponding prefractal flat surface S{KSn). In §2.1 we saw how to 
construct a flat surface from a rational billiard table. In the case of the equilateral 
triangle billiard table r2(A) = ^{KSq), there are 2 • lcm{3, 3, 3} = 6 copies of f2(A) 
used in the construction of the associated flat surface 5(A); see Example 2.6 and 
the associated Figure 2. In the case of the prefractal billiard table il{KSn), only 
six copies of r2(i^5'„) are needed in the construction of the associated flat surface 
S{KSn), for every n > 0; see Figure 19. (We refer to [LapNiel, LapNie2, LapNie3] 
for further discussion of the topics in the present subsection.) 

The vertices of il{KSn) correspond to conic singularities of the flat surface. How- 
ever, only certain singularities are removable. The vertices with angles measuring 
^ (measured from the interior), constitute removable singularities of the flat sur- 
face. That is, the geodesic flow can be appropriately defined at these points. The 
vertices with angles measuring ^ constitute nonremovable singularities. Hence, it 
is possible to define reflection at certain vertices of the prefractal billiard il{KSn), 
but impossible to define at others. Moreover, defining reflection at acute corners 
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of fl{KSn) in this way is independent of n. That is, how a biUiard ball reflects in 
an acute corner — as dictated by the corresponding geodesic on the associated flat 
surface — never changes from one prefractal billiard table to the next. 

Such insight is clearly helpful in further understanding the behavior of a billiard 
ball on the Koch snowflake fractal billiard fl{KS), but we must be careful not to 
extrapolate more than is possible from this observation. Knowing that we can 
determine an orbit of a prefractal billiard n{KSn) by unfolding the orbit of n{KSo) 
in n{KSn), we are inclined to allow orbits of ^{KSq) that make collisions with 
corners. However, a priori, we cannot conclude that such orbits do not unfold to 
form saddle connections in n{KSn) connecting two nonremovable singularities. In 
the event an orbit of rt{KSm) intersects the boundary KSm solely in 

acute corners, then such an orbit is an element of a sequence of compatible orbits 
with ^j{x°,e°) = for every j > m. 

4.3. The T-fractal prefractal billiard. We refer to §3.3 for a discussion of the 
T-fractal ^ and of its prefractal approximations £^n, for n = 0,1,2,...; see, in 
particular. Figure 12. Recall that the base of has a length of two units. The 
prefractal billiard ri(5^) can be tiled by the unit square Q; see Figure 13. In 
general, for every n > 0, ^{^n) can be tiled by the square ^Q. As such, and since 
Q obviously tiles the plane, we can apply Theorems 4.7 and 4.9. 

Much like the case of the prefractal Koch snowflake billiard il,{KSn), we are in- 
terested in forming sequences of compatible orbits of prefractal billiards exhibiting 
particular properties. The results in this subsection appear here for the first time 
and will be further discussed in [LapNie6] . We recognize that if a periodic orbit has 
an initial condition (ccq,^})), then there may exist a compatible orbit ffN{x%,9%) 
that forms a saddle connection if ajg has a finite binary expansion. This is not to 
suggest that ffN{x%, ^^v) n^'^^st form a saddle connection. However, Xq having an in- 
finite binary expansion (and no equivalent finite binary expansion^^) and ^0(2^0, ^0) 
being a periodic orbit guarantees that {^„(a;^, ^'')}^o ^ sequence of compatible 
periodic orbits. 

Theorem 4.23. Let {xq,9q) be an initial condition of an orbit i^o(a;o,^o) ofil{Q). 
If every element of the footprint J^o{xq,9q) only has an infinite binary expansion 
{and no equivalent finite binary expansion), then the sequence of compatible orbits 
{^nix"^,0^)}^^o {where {x'^,e^) = {x'^,e^)) of the prefractal billiards n{X) is a 
nonconstant {or more precisely, is not an eventually constant) sequence of compat- 
ible of periodic orbits. 

Example 4.24. Let a;g = | and 6'g = f . Then, {^„(a;°, f )}~=o is a nonconstant 
sequence of compatible periodic orbits; see Figure 20. 

Theorem 4.25. Let Xq = with fc, t e N, t and 3 relatively prime, k ^ and 
< f < 3*^. Further, let m e M. If for every p, q,r,s G Z, r,s > 0, we have that 

nr-sok 

(14) Z , 

then the line y = m{x — Xq) does not contain any point of the form 
a, b,c,dG Z, with c,d> 0. 



^^Recall that a value may have two equivalent binary expansions. For example, i = 0.1 = 0.01. 
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Figure 20. A sequence of compatible periodic orbits of i}{^o), 
Sl{^i) and fl{^2), respectively. 




II 12 21 ll 1: 12 21 n 11 12 21 n 



Figure 21. The flat surfaces <S(,5?)), Si-^^i) and associated 
with the T- fractal prefractal approximations and 
respectively. 

Note that the condition (14) above is automatically satisfied if the slope m is 
irrational. 

Theorem 4.26. Let {xQ,yQ) = (p-,0), with k,t G N, t and 3 relatively prime, 
k ^0 andO <t <3''. If 

27 

(2a+l)^' 

with a, /3, 7 € N, a,(3,j > 0, then, for every p,q,r,s G Z with r,s > 0, the point 
^) does not lie on the line y = m(x — Xq). 

4.3.1. The corresponding prefractal flat surface S{^n)- For every n > 0, the in- 
terior angles of ,%i are j and To form the associated flat surface S{,%i), we 
appropriately identify four copies of ^{3^n)', see Figure 21 for a depiction of the first 
three flat surfaces. 

Then, every point of <S(,5^) associated with a vertex of r2(,5^) measuring ^ con- 
stitutes a removable singularity of iS(=^). Similarly, every point of iS(^) associated 
with a vertex of n(^) of interior angle measuring ^ constitutes a nonremovable 
singularity of S{.%i). Therefore, not every vertex of 0,{^n) will present a problem 
for the billiard flow. 

Consider an orbit of fl{Q), where the orbit has basepoints corresponding to ver- 
tices of Q, the unit square. Since such vertices correspond to removable singularities 
in the corresponding flat surface (this being the flat torus, see §2.1), we see that the 
same orbit reflected- unfolded in the billiard O(^) (if one flrst scales the billiard 
f2((5) and the orbit contained therein by ^rr, see §2.2) can potentially intersect ver- 
tices of that are associated with nonremovable singularities in the corresponding 
flat surface. 
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4.4. A prefractal self-similar Sierpinski carpet billiard. Let Sa be a self- 
similar Sierpinski carpet, as defined in Definition 3.8, and let us denote its natural 
prefractal approximations by Sa,i for z = 0, 1, 2, ... (as in §3.4). The corresponding 
billiard is then denoted by ^l{Sa)- In this subsection, we examine the behavior of 
the billiard flow on the rational polygonal billiard given by the prefractal approxi- 
mations fl{Sa,i)^'^ In the event a billiard ball collides with a corner of a peripheral 
square, we must terminate the flow and such a trajectory is then called singular. In 
addition to being singular, such a trajectory will form a saddle connection (see the 
beginning of §2 for a discussion of closed billiard orbits that form saddle connec- 
tions). As we have discussed, an examination of the corresponding flat surface may 
prove useful in determining whether or not a billiard ball can reflect in a vertex. 

Definition 4.27 (Obstacle of fl{D)). Let fl{D) he a polygonal billiard. Then 

fl(D) can be modified by placing in its interior a piecewise smooth segment that 
inhibits the billiard flow and causes a billiard ball to reflect. Such a segment is 
called an obstacle of Q{D). 

Clearly, each prefractal billiard Cl{Sa,i) can be interpreted as a square billiard 
with obstacles. 

Notation 4.28. Due to the fact that Theorem 3.13 refers to the slope of a nontrivial 
line segment and we make heavy use of this theorem, we will denote the initial 
condition {x^, e^) of an orbit of 0(5a,„) by (a;° , ), where = tan(6»°). 

Definition 4.29 (An orbit of the cell Ck^a'' of ^{Sa,k))- Consider the boundary 

of a cell Cj, (jfc of n{Sa.k) as a barrier. "'^^ Then an orbit with an initial condition 
contained in the cell is called an orbit of the cell Cf^ ^'' of n{Sa,k)- 

Remark 4.30. So as to be clear, the boundary of the cell does not form an obstacle 
to the billiard flow, as deflned in Definition 4.27. Rather, we are treating the 
cell Cj, „fc as a billiard table in its own right, embedded in the larger prefractal 
approximation il{Sa.k)- 

Recall from §3.4 that a self-similar Sierpinski carpet Sa is the unique fixed point 
attractor of a suitably chosen iterated function system consisting of simi- 

larity contractions. In light of this, an orbit of a cell C^, ,jfc of n{Sa,k) is the image of 
an orbit &o{xq, Uq) of the unit-square billiard fl{So) under the action of a composi- 
tion of contraction mappings (pmk o • • • o 4'mi j with 1 < rrii < — 1 and I < i < k, 

2_i 

determined from the iterated function system {(t>j}j^i of which Sa is the unique 
flxed point attractor. 

Lemma 4.31. Consider a self-similar Sierpinski carpet Sa- Let fc > and Sa.k 
be a prefractal approximation of Sa- If a £ Ba,^''^ then the line segment beginning 
at a midpoint of a cell C^^a'' of Sa,k is a nontrivial line segment {in the sense of 
Definition 3.12). Moreover, such a segment avoids the boundary of the peripheral 
squares of Sa with side-length a"™, m > fc -|- 1. 



We note that the results in this subsection appear here for the first time and will be further 

discussed in [CheNie|. 

^^Hcre, Cf, is a cell of the feth prefractal approximation Sa,kt 3* given in Definition 3.10 
with all numbers aj equal to a. 

^"Recall from Notation 3.14 that Ba is the set of slopes given by Equation (13). 
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Figure 22. Interpreting the flat surface S{Sa,n) as a flat torus 
with obstacles. 



The statement in Lemma 4.31 assorts that a segment beginning at a midpoint 
of a cell with slope a & Ba will be a nontrivial line segment in Sa- In addition to 
this, any line segment contained in M? that contains a nontrivial line segment of Sa 
must necessarily avoid the peripheral squares in a tiling of by Sa- Otherwise, 
there exists k > 1 such that scaling the line segment in and the tiling of by 
a"'' results in a segment contained in the nontrivial line segment which intersects 
peripheral squares of Sa- This is a contradiction of the fact that the segment 
beginning at (2~^, 0) with slope a e is a nontrivial line segment of Sa- We then 
deduce the following result. 

Theorem 4.32. Consider a self-similar Sierpinski carpet Sa- Let k > Q and Sa,k be 
a prefractal approximation of Sa- Furthermore, let a € Ba and = (p(2a'^)~-^, 0) 
with p < a positive, odd integer. If i7'k{x^,a^) is an orbit of ^l{Sa,k), then 
the initial condition (a;^,a^) determines a sequence of compatible periodic orbits 
{ffn{Xnj <^^)}^=k of the prefractal approximations fl{Sa,n)- 

As one may suspect, there exists N > k > such that a sequence of compati- 
ble orbits {^„ (.t" , a'^)}^^pf is a constant sequence of compatible orbits. Moreover, 
— ^'n' for every n > N- This is not any different from the case of a constant 
sequence of compatible orbits of prefractal billiards Cl{KSn), as discussed in The- 
orem 4.21 and Example 4.22. However, in the context of a self-similar Sierpinski 
carpet billiard table, every sequence of compatible orbits we will examine will be 
a sequence for which there exists N >0 such that {^n{XnT oP)}'^^=n ^ constant 
sequence of compatible orbits. 

4.4.1. The corresponding prefractal flat surface S{Sa,i)- In much the same way the 
billiard Q{Sa,i) can be interpreted as a square billiard with obstacles, the corre- 
sponding flat surface can be interpreted as a "torus with obstacles"; see Figure 
22. 

In light of the fact that S{Sa,n) can be interpreted as a torus with obstacles 

and the presence of a dynamical equivalence between the billiard flow and the 
geodesic flow on the corresponding flat surface (see §2.2), we see that reflection in 
the vertices with angles measuring ^ (relative to the interior) can be deflned. More 
specifically, the geodesic flow can be defined at points corresponding to vertices 
with angles measuring 7r/2, because these points constitute removable singularities 
of the geodesic flow. 
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This fact is crucial in determining orbits of 51 (5a) for which the slope a. is an 
element of Aa and not Ba (see Notation 3.14), and the orbit avoids all peripheral 
squares of ^(S'q). While the reader may say that this contradicts Theorem 3.13 
(and they would be right), in [CheNie] a more precise formulation of Theorem 3.13 
is given that clarifies which slopes are permissible and which ones are not. That 
is, if a° e Aa, it may be possible for an orbit ^„(x",a") to begin at the origin 
and avoid the peripheral squares of each billiard ^{St^m), for every m > n and 
21 + 1 = & < a, for some / > 1. We do not give here an explicit reformulation of 
Theorem 3.13, but Example 5.12 in §5.4 exhibits a situation showing that Theorem 
3.13 is not stated sufficiently precisely. 

5. Fractal billiards 

The theme that will tie together all of the examples in §4 is that suitable limits 
of sequences of compatible orbits may constitute billiard orbits of each respective 
fractal billiard table. We have shown that in the case of il.{KS), J7(,5') and r2(5a), 
we can determine a sequence of compatible periodic orbits. We will see that in 
each case of a fractal billiard, under certain conditions, a sequence of compatible 
periodic orbits (or a proper subset of points from each footprint J>i(x",(?")) will 
converge to a set which can be thought of as a true orbit of a fractal billiard table 
(or such a sequence will yield a subsequence of basepoints converging to what we 
are calling an elusive point in [LapNie2, LapNie3]). 

5.1. A general framework for il.{KS), and rt{Sa)- We restrict our at- 

tention to the family of fractal billiard tables fl{F) where F is a fractal approxi- 
mated by a suitable sequence of rational polygons {F„}^q, with each tiled by 
Dn = CnDf) for suitably chosen c„ € (0, 1] and Dq a polygon that tiles the plane. 
Specifically, we are interested in developing a general framework for dealing with a 
fractal billiard table 0(F) that is constructed in a way which is similar to that of 
Vl{KS), Vt{.^) and Vt{Sa). 

Before we begin our discussion of the fractal billiard tables ^1{KS), and 
Q.{Sa), we define certain terms. The following definitions were initially motivated 
by the work in [LapNieS], but later generalized for this paper in order to account 
for a larger class of fractal billiard tables. (From now on, we assume that ^{F) 
is a fractal billiard table with prefractal billiard approximations {f^(-Fn)}^o ^ 
described just above.) 

Definition 5.1 (A corner). Let z € F. If there exists n > such that z € Fn and 

2; is a vertex of Fn- then z is called a corner of F. 

Definition 5.2 (A Cantor point). Let z G F he such that z is not a corner of 
F. If there exists N > such that for every n > N, z £ Fn and every connected 
neighborhood of z contained in F„ becomes totally disconnected when intersected 
with F, then z is called a Cantor point of F. 

In the Koch snowflake KS, every Cantor point is a smooth point of infinitely 
many prefractals KSn approximating KS. That is, if 2: is a Cantor point in KS, 

then there exists iV > such that for every n > N, there exists a well-defined 
tangent at z e KSn-^'^ We deduce from this that the law of reflection holds at 



Here and Definition 5.3 below, 2 is viewed as a point of the smooth subarc of F„ to which it 
belongs. 
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z e KSnt for every n > N. Moreover, since the billiard ball reflects at z G KSn at 
the same angle for every n > N, we deduce that the tangent at z is the same for 
each KSn, n> N. This observation then prompts us to generalize the definition of 
a Cantor point in order to account (for example) for points of the T-fractal which 
are not Cantor points, but are points for which a well-defined tangent can be found 
in infinitely many prefractal approximations. 

Definition 5.3 (Smooth fractal point). Let z & F and iV > be such that 2; G F„ 
for every n> N. If there exists a well-defined tangent z & Fn for every n> N, 
then z is called a smooth fractal point. 

Definition 5.4 (An elusive point). Let z G F. li z ^ U^o ^^^'^ ^ called an 

elusive point of F. 

Consider a piecewise linear path in ^{F), such that every linear segment of the 
path is joined at the endpoint of another segment with the coincidental endpoints 
intersecting the boundary _F at a smooth fractal point of F (in the sense of Definition 
5.3). In the following definition, we define a particular type of piecewise linear curve 
in a fractal billiard Q.{F). 

Definition 5.5 (A nontrivial path). Suppose that there exists a piecewise linear 

curve in Vi{F) as described immediately above. If at each point z for which the 
piecewise linear path intersects the boundary F, the angle formed by the first 
segment is equal to the angle formed by the second segment, relative to the side of 
Fn on which z lies,^^ then the piecewise linear path is called a nontrivial path of 

n{F). 

Remark 5.6. In [LapNieS], a nontrivial path was called a nontrivial polygonal 
path. The change in name is purely based on aesthetics. 

Definition 5.7 (A Cantor orbit). Suppose ^„(a;°,6'°) is an orbit of Q.{Fm), for 
some > 0, such that every point of the footprint J^n{x%,9%) corresponds to a 
smooth fractal point of F. This then readily implies that ^„(a;° , 0°) is the same as 
^n{-i:%,0%) for every n > N.^"^ Then ^„(a;°,6»°) is called a Cantor orbit of 0(F) 

and is denoted by ff{x^,9^). 

If O(i^) is a fractal billiard table, then it may or may not be possible to construct 
Cantor orbits or nontrivial paths of fl{F). We will next discuss three examples of 
fractal billiard tables with different dynamical properties that lend themselves well 
(or not) to determining well-defined billiard orbits. 

Remark 5.8. We note that applying Definitions 5.1, 5.2 and 5.4 to n{KS) and 
the sequence of rational polygon prefractal approximations 0(A'5'„) which we have 
discussed in §3.2 yields exactly the sets of points we are considering as corners, 
Cantor points and elusive points of fl{KS), respectively. Moreover, applying Defi- 
nitions 5.1, 5.3 and 5.4 to O(^) and the prefractal approximations Ct{^n) which we 
discussed in §3.3 yields exactly the sets of points that we arc considering as corners, 
smooth fractal points and elusive points of Finally, applying Definitions 5.1 

and 5.3 to il{Sa) and the prefractal approximations il{Sa^n) which we discussed 

^^Recall from Definition 5.3 that a smooth fractal point 2 of -F is a point of infinitely many 
prefractal approximations Fm- Hence, there is a least nonnegative integer n such that z G Fn- 

^^In other words, {fin{Xn, 6^)}^—isr is a constant sequence of compatible orbits, where a se- 
quence of compatible orbits was defined in Definition 4.5. 
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in §3.4 yields exactly the set of points we are considering as corners and smooth 
fractal points of Cl{Sa)- 

5.2. The Koch snowflake fractal billiard. As we have noted before at the end 

of §3.2, for each n > 0, KSn O KS can be realized as the union of 3 • 4" self-similar 
ternary Cantor sets, each spanning a distance of Within each Cantor set, we 
find Cantor points and corners of the Koch snowflake. 

We begin our discussion of orbits of fl(KS) by examining the limiting behavior 
of a particular sequence of compatible orbits with the initial condition {x%, J), 
where is a Cantor point of KS (i.e., x% is a point of KSn with a well-defined 
tangent in KSn for every n > N). For the sake of simplicity, we let TV = and 
= i be on the base of the equilateral triangle KSq (recall that we are assuming 
that the left corner of KSo is at the origin and the length of each side is one 
unit). Then, ^o(a"o-f) is an orbit that remains fixed as one constructs n{KSi) 
from Q{KSo). More correctly, {^„(a;° , f )}J^o ^ sequence of compatible orbits 
with J^n{Xn, f ) = ^o{xq, f ) for every n > (that is, with the same footprint in 
each prefractal approximation). 

In general, if (cc^, ^) is an initial condition of an orbit of ^{KSn) and 
Cantor point, then the sequence of compatible orbits is such that for every n> N, 
the footprints J>i(a;^, ^) and J^n{x%, |) are the same. 

Theorem 5.9. If x € KS is a Cantor point, then there exists a well-defined orbit 
of fl{KS) with an initial condition (x,-n7(|-)), where the angle a7(^) is determined 
with respect to the side on which x lies in a prefractal approximation Q,{KSn)- 

There are many more properties of {^„(a;° , f )}^jv which we could discuss here. 
These properties largely rely on the nature of the ternary representation of x^ , and 
are elaborated upon in [LapNie2, LapNie3]. We now proceed to illustrate how we 
can connect two elusive points of Sl{KS). Such a result has already been presented 
in greater detail in [LapNie3], so we will be brief. In §5.3, we will show that an 
identical construction holds for the billiard table f2(^). 

Recall from Example 4.20 that we were able to construct a sequence of compatible 
periodic hybrid orbits. From such a sequence we can derive a sequence of basepoints 
that is converging to an elusive point of Q{KS). The latter sequence of basepoints 
constitutes the vertices of a nontrivial path; see Figure 23. One may consider 
a direction 7q that is the reflection of through the normal at Xq. Then, the 
resulting sequence of compatible periodic hybrid orbits {^n{Xn^lo)}^=N yields a 
sequence of basepoints converging to another elusive point. Again, such a sequence 
of basepoints constitute the vertices of a nontrivial path of rt{KS); see Figure 24. 
Together, these two nontrivial paths constitute a single nontrivial path connecting 
two elusive points of fl{KS). 

In conjunction with Theorem 4.19, we can determine countably infinitely many 
initial conditions {x^, 0^), each of which determines a sequence of compatible peri- 
odic hybrid orbits yielding a sequence of basepoints converging to an elusive point 
of n{KS). 

5.3. The T-fractal billiard. The results in this subsection appear here for the 
first time and will be further discussed in [LapNie6]. We begin our discussion of 
the billiard by recalling (and referring the reader back to) Example 4.24 from 
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Figure 23. A nontrivial path of the Koch snowflake fractal bil- 
liard table Q{KS) beginning at a; = ^. 




Figure 24. Two nontrivial paths connecting two elusive points 
of n(KS). (As is explained in the text, these two paths can be 
concatenated to obtain a single nontrivial path connecting the two 
elusive points.) In the first figure, we only show the relevant por- 
tions of the Koch snowflake. In the second figure, we magnify the 
regions containing the nontrivial paths so as to highlight the fact 
that such paths are converging to elusive points. Actually, there is 
an obvious geometric similarity one can take advantage of in order 
to produce more segments of the nontrivial path. 




Figure 25. Two nontrivial paths connecting two elusive points of 



§4.3. The sequence of compatible periodic orbits provided by Example 4.24 gives 
rise to a nontrivial path that connects | with an elusive point of Further- 
more, considering the sequence of compatible periodic orbits {^„(|, we 
determine another nontrivial path that connects | with another elusive point of 
ri(c^J'); see Figure 25. This behavior is analogous to the one which we observed for 
the Koch snowflake billiard in §5.2. 

As was the case with ^1{KS), we can analogously build upon Theorems 4.25 and 
4.26 in order to determine a sequence of basepoints converging to an elusive point. 
That is. Theorems 4.25 and 4.26 guide our search for a sequence of compatible 
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periodic orbits which yields a sequence of basepoints converging to an elusive point 

Theorem 5.10. Let {^n{Xn,d^)}^^N ^ sequence of compatible orbits. Then, 
there are countably infinitely many directions and countably infinitely many points 
from which to start in order to produce two nontrivial paths that will connect two 

elusive points o/0(,5^). 

5.4. A self-similar Sierpinski carpet billiard. In [Du-CaTy], nontrivial line 
segments of Sierpinski carpets are constructed. Building on the main results of 
[Du-CaTy], the second author and Joe P. Chen have been able to construct a 
family of periodic orbits of a self-similar Sierpinski carpet, in the sense of [LapNie2, 
LapNie3] recalled in Definition 5.7.^"' Such orbits constitute Cantor orbits of the 
self-similar Sierpinski carpet. As of yet, we have not attempted to construct a 
nontrivial path of a Sierpinski carpet. 

In light of Theorem 4.32, we say that the trivial limit of a constant sequence of 
compatible periodic orbits constitutes a periodic orbit of a self-similar Sierpinski 
carpet billiard n{Sa)- In the event an orbit has an initial direction aj], we may 
still be able to determine a constant sequence of compatible periodic orbits. The 
trivial limit of such a sequence then constitutes a periodic orbit of Q{Sa)- Using the 
fact that reflection can be defined in the vertices with interior angles measuring J, 
we can state the following result. (Recall from §4.4 that Sa,n is the nth prefractal 
approximation of Sa-) 

Theorem 5.11. Recall from Notation 4-28 that if is the initial direction of a 
billiard orbit, then a = tan^. Let x° = (0,0), a € Q and let i)'{x^,a) be an orbit 
offl{So)- If &{x'^,a), as an orbit of il{Sa,i), avoids the middle peripheral square, 
then the initial condition {x°,a) will determine an orbit ofO,{Sa)- Specifically, the 
path traversed by the orbit ff{x'^,a) of fl{Sa,i) is exactly the path traversed by the 
orbit ofn{Sa) determined by (x'',a). 

Example 5.12. Let x° = (0,0), a = 2/3 G Slope(S'5). Consider an orbit of 
^{87,2) with an initial condition {x^,a); see Figure 26. We see that the orbit 

avoids the peripheral square of Q{St_i). By Theorem 5.11, the initial condition 
{x'^,a) determines an orbit of Q{Sr). The path traversed by the orbit of fl{Sr) is 
exactly the path traversed by the orbit ^{x^, a). 

6. Concluding remarks 

It is clear from the preceding sections that much work remains to be developed 
in order to determine a well-defined phase space {^{F) x S^)/ ^ for the yet to be 
defined fractal billiard flow. We have discussed several examples of what clearly 
constitute periodic orbits of Cl{KS) and 0,{Sa)- Furthermore, for both i}{KS) and 
we were able to connect two elusive points of each billiard table via suitably 
chosen nontrivial paths. These nontrivial paths were determined from suitably 
chosen sequences of compatible periodic orbits. 



The results described in this subsection appear here for the first time and will be further 
discussed in [CheNie]. 
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Figure 26. An orbit with an initial condition beginning at (0,0) 
and with an initial direction constituting a slope of a = 2/3 e Ay, 
where is defined as in Notation 3.14. While it would appear 
that this orbit intersects corners of peripheral squares, it in fact 
remains away from all peripheral squares. The same is true for 
finer approximations. 

Question 6.1. Let F be either KS or ^. Suppose that two suitably chosen non- 
trivial paths converge to two distinct elusive points ofil{F). For each of the two 
elusive points, is it possible to determine another nontrivial path converging to a 
different elusive point? 

If we can answer Question 6.1 in the affirmative (or answer it in the affirmative 
under specific conditions), will this help us gain insight into how to determine 
a well-defined phase space for the billiard flow on n{F)7 An alternate approach, 
discussed in the concluding remarks of [LapNie3] , entails determining a well-defined 
fractal flat surface. 

Question 6.2. Regarding a self-similar Sierpinski carpet billiard Cl{Sa), we have 
determined a countable set of points from which a periodic billiard orbit can begin. 
Can we show that the set of points from which a periodic orbit can begin is in fact 
uncountable and, furthermore, a set of full {Lebesgue) measure in the base of the 
unit square So ? 

Question 6.3. In analogy with the prefractal billiard and associated flat surface, 
can a thorough understanding of the geodesic flow on the limiting {and still to be 
mathematically defined) ^fractal flat surface^ '5(F) aid us in determining a well- 
defined billiard flow on fl{F) ? 

The work in progress in [LapNic4] draws upon the work of Gabriela Weitze- 
Schmithiisen [We-Sc] and attempts to answer Question 6.3 from an algebraic per- 
spective. 

Approaching the problem of determining a well-defined billiard flow on a fractal 
billiard table from many different angles may prove useful. The theories of flat 
surfaces and rational billiards are intimately tied together and more deeply under- 
stood by knowing the structure of what is called the Veech group (this being the 
group studied in [We-Sc, Ve3, Ve4, Vo]). In short, the Veech group of a flat surface 
S{D) determined from a rational polygon D is the stabilizer of S{D). 

Question 6.4. Let fl,{F) be a fractal billiard table, with F being approximated 

by a suitably chosen sequence of rational polygons {FnlJ^Q. Js it then possible to 
construct a Veech group for ^1{F) {or rather, of S{F)), presumably in terms of the 
Veech groups for the prefractal approximations f2(F„) (or rather, of the associated 
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flat surfaces Will the knowledge of such a group aid us in determining a 

well-defined billiard flow on 0(F) ? 
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